Entanglement Localization and Optimal Measurement 
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The entanglement can be localized between two noncomplementary parts of a many-body system 
by performing measurements on the rest of the system. This localized entanglement (LE) depends 
on the chosen basis set of measurement (BSM). We derive here a generic optimality condition for 
the LE, which, besides being helpful in studying tripartite systems in pure states, can also be of 
use in studying mixed states of general bipartite systems. We further discuss a canonical way 
of localizing entanglement, where the BSM is not chosen arbitrarily, but is fully determined by 
the properties of the system. The LE obtained in this way, we call the localized entanglement 
by canonical measurement (LECM), is not only operationally meaningful and easy to calculate in 
practice (without needing any demanding optimization procedure) , it provides a nice way to define 
the entanglement length in many-body systems. For spin-1/2 systems, the LECM is shown to be 
optimal in some important cases. At the end, some numerical results are presented for ji — ]2 spin 
model to demonstrate how the LECM behaves. 

PACS numbers: 03.67.Mn, 64.70.Tg, 03.67.-a 



Besides its importance in interpreting and understand- 
ing quantum mechanics, the entanglement has gained im- 
mense interest in recent times as it has the potential to 
play a significant role in modern technology. In addi- 
tion, it has become an important tool to study quantum 
many-body systems Some of the very useful mea- 
sures used here for studying quantum non-local nature of 
a system are pairwise entanglement [U, 0] j local entropy 
[if, localizable entanglement @ and negativity @. Here 
we study entanglement localization which is important 
for two reasons -from a practical point of view, it can be 
a useful method of producing entangled pairs (especially 
from three-body systems) and secondly, we get an alter- 
native theoretical way for studying quantum many-body 
systems. Here it may be stressed that, almost all the 
measures which try to quantify mutual quantum behav- 
ior between two disjoint parts of a many-body system are 
cither difficult to calculate (often needing optimization 
procedure) or they do not have any operational mean- 
ing. The measure we present here (LECM) is not only 
easy to calculate, it also has some operational meaning. 

Let S\ and £2 be any two noncomplementary parts of a 
total system U . The rest of the system is called the envi- 
ronment (E), which generally consists of many sites (Fig. 
[LJl) . A measurement on E by some basis set would result 
in S (= Si + S2) assuming different pure states with ap- 
propriate probabilities. Unlike in the case of localizable 
entanglement Q where only local measurements on the 
individual sites of E are allowed, we allow all possible 
measurements (including the joint measurements on the 
sites) in our localization process. It may be noted here 
that, all the measurements in this work arc considered 
to be non-selective projective-type. In the next part, we 
derive a simple but generic optimality condition for the 
LE, which will help us find optimal (more generally, sta- 
tionary) solutions and check whether a given solution is 
optimal. Studying a general bipartite system in mixed 
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(b) Chain 

FIG. 1. (a) Total system or Universe (U): Si and S2 are two 
noncomplementary parts (they need not be identical) that 
together form the system S. Rest of U is the environment E. 
(b) Chain: the symmetrically placed (about dotted C2 axis) 
filled circles represent two parts/sites under study. 



state is notoriously difficult. There can be innumerable 
ways of decomposing a mixed state, where each decompo- 
sition corresponds to an average entanglement (entropy). 
The maximum and the minimum possible values of the 
average entanglement are termed as entanglement of as- 
sistance (EoA) [3] and entanglement of formation (EoF) 
[H respectively. The optimality condition derived here 
may be of use in finding them. In this regard, a brief 
discussion is given after arriving at the condition. 

The optimality condition.- When expressed in the 
product basis states of E and S (Fig. OJi), the 
given wave function (that we study) becomes, = 

£5=5 C t M) E \^) S - Here |0 B s (|0) s s) arc some or- 
thonormal basis vectors of the state space of E (S) with 
dimensionality De (Ds)- The state can also be written 
as 



z — 1 



(i) 



with Pi = v<: : (;., (■;.,■ and = Y,%%\<l>i) s - 

Here the summation runs over nonzero Pi's, numbering 
D (< De)- In general, states \£) s s, are not orthonormal. 
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The operational interpretation of the later expression of 
the state |^) is that, if we perform measurement on E 
by the basis set {£ E }, the state will collapse and we will 
get S in different pure states \£,i) s s with corresponding 
probabilities p^s. 

If Si be the entropy of , then the average entropy 
(entanglement) localized between Si and Si would be, 

S{^} = £%A. (2) 

As both p^s and Si's depend on the choice of the BSM, 
the average entropy (or LE) S will also depend on the 
choice of the BSM. We need to derive a condition for the 
choice of BSM ({^ s }) which optimizes S. 

We first note that, any general basis set can be ob- 
tained from an initial basis set {£, E } by application of 
a series of elementary transformations (ETs). Here an 
ET is a small-angle orthonormal transformation between 
any two initial basis states keeping others unchanged. We 
now derive first order change in S due to an ET. If 
and \£,j) E be any two initial basis states, then the two 
new basis states obtained by an ET would be, 

Wi) E = \& E + z\*Li) E and \i') E = \Q E e\d t ) E . (3) 

Here e is the small angle (a parameter) whose higher 
order terms can be neglected. Due to change in these 
basis states, corresponding probabilities and states of the 
S would also change (see eq. [1]). We need to relate these 
new probabilities and states with the old ones. 

At this stage it is advantageous to express all the 
probabilities as the diagonal elements of a density op- 
erator (matrix), which is, in our case, the reduced den- 
sity matrix (RDM) of E (p E ). The elements of the 
RDM are given by p E = V" . (\ ,C; ,. Using this 
RDM, probability corresponding to a state \£) E would 
be p = E (£\p E \£,) E ■ This allows us to write the new 
probabilities as (using eq. [3]), 

Pi = Pi + e hj and p rj = pj - efcjj, (4) 

with % = k Tl = E (^\p E \^) E + E (^\P E \^) E - Let us 
first consider the case when none of the pi and pj is zero. 
Now if |£-) 5 and be the new states of the S, then, 
in the new scenario, the state \^f) can be rewritten as, 

i*> = M& E \& S + ftj\® E \® s + ■■■ (5) 

Here we focus only on i-th and j-th states, as other states 
arc unchanged. Now using cqns. [3] and HI in the above 
expression and then comparing the terms associated with 
the initial basis states \£,i) E and \£j) E from the two dif- 
ferent expressions of (in cqs. [T]and[5]), we get the 
following solutions for the new states of S: 

lO s He t ) 5 + e K-i&) 5 + M^) 5 ) (6) 
iO s H&> s -eM&> s +M&> s ) (7) 



Here = — \kijp { 1 and bij = p] Pi in eq. [H In- 
terchanging the indices i and j we get similar terms in 
eq.El 

Now let Q = {Qi m } and R = {Ri m } be the matri- 
ces representing respectively the states \£i) s and \£j) s 
in some product basis states of the parts Si and 5*2. 
In terms of these matrices, the RDMs for Si would be 
p Sl (£_i) = QQ^ and p Sl (£j) = RR f while S is respec- 
tively in \&) s and \^) s . Similarly, the RDMs for Si 
corresponding to the new states, given in eqs. [5] and 
would be, 

P Sl (£) = P Sl (fi) + e {2a ijP Sl (6) + 2b ij A ij ) (8) 
P Sl - P Sl ~ e (2«ii/ J &) + (9) 

Here A !3 - = ^(QR^ + RQ^), a Hermitian matrix. Let us 
denote here the changes in the RDMs in cqs. [5] and |H] as 
epi % (ij) and —ep^^i) respectively. It is worth mention- 
ing that, as trace (Tr) of any RDM is 1, we have 

Trpf (ij) = Trpf (jj) = 0. (10) 

We now use the relation p Sl (£■) log2 p Sl (CD = 
log 2 ^(&) + epf^ij) + epp{ij) log 2 p Sl (&) 
for obtaining entropy corresponding to the new state 
\d) S ■ This operator relation is not ill-defined due to last 
term as both p Sl (£i) and p^iij) go to zero simultane- 
ously (this can be understood by singular value decom- 
position of the matrix Q). Now tracing over both sides of 
this relation and a similar relation for the j-th state, we 
respectively get the following entropies for the new states 
of S, 

Si = Si - eTr pf 1 (ij) log 2 and (11) 

5J = S 3 + eTr pf 1 (ji) log 2 p s > (fc ) . (12) 

Here we used eq. [TU] to get these relations. Let us now 
denote the changes in entropies in eqs. [11] and [T2] as — eSy 
and eSji respectively. 

Now using these new entropies along with the new 
probabilities in eqn. 21 we get the new average entropy: 

S' = Y^ =1 P'iSl = S + eSi, (13) 

where, Si = fc^S; - p<5y - kjiSj + PjS^. 

Before we set the optimality condition, we now check 
the cases when both are or one of pi and pj is zero. When 
Pi = pj = 0, then kij = kji = 0. Therefore, p[ = p'j = 
(see eq. H]). Which implies that Si is zero. On the 
other hand, when pi ^= and pj = 0, we gave again 
kij = kji = 0. From eq. [4] we have p\ = pi and p'j = 0. 
Now it is clear from eq. [5] (with second term being zero) 
that, this type of ETs are not allowed (within the first 
order calculation). 

So, the desired optimality condition is Si — or, 

hj Si — piSlj = kjiSj — PjSji, (14) 
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for all i and j for which corresponding probabilities are 
nonzero. The second order change in S due to different 
ETs can also be derived but they can not confirm actual 
character of an optimum flpj ] . 

For the derivation of the above optimality condition, 
we have assumed the existence of a fixed environment 
(E). Therefore, the condition will be helpful when we 
have a definite tripartite system and we want to local- 
ize entanglement between two parts in an optimal way 
by performing measurement on the third part. Now the 
question is whether it also can be of any use for calculat- 
ing EoA and EoF of a mixed state. We here note that, 
for a bipartite system in a mixed state, it is always possi- 
ble to construct a pure state by augmenting the bipartite 
system with an ancilla in such a way that the RDM of 
the system becomes the given mixed state. By perform- 
ing all possible measurements on all possible ancillas, we 
get all possible decompositions of the mixed state. By a 



theorem of Hughston- Jozsa-Wootters [11( , we can relate 
each decomposition to an M rX k matrix with k orthonor- 
mal column vectors. Here the r and the k are the number 
of terms in the decomposition and the rank of the mixed 
state respectively. This implies that, for a particular an- 
cilla, we can express the optimality condition in terms of 
a matrix M rX k- Now as in principle number of terms in a 
decomposition can be anything, it may appear that the 
optimality condition obtained for a fixed ancilla would 
be of no use for the said purpose. Fortunately, at least 
in the case of EoF, it is seen that, consideration of a 
few number of terms in the decomposition is enough for 
Hence we hope that, optimality con- 
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extremization 

dition given here may also be useful in calculating EoA 
and EoF of a mixed state. 

The LECM.- Our canonical way of localizing entangle- 
ment between S\ and S2 is to take eigenstates of p E (the 
RDM of E) as the BSM and perform measurement on 
E. We will find now the expression for the entanglement 
localized in this way (which we call LECM). 

Expression of the LECM can easily be obtained from 
Schmidt decomposition (SD) of the state under study. 
The SD of the state into the product states of E and 
S is given by, 



I*) 



(15) 



Here, \£.) s (\£.) E ) is the i-th eigenstate of the RDM of 
S (E) corresponding to the eigenvalue A^. The index i 
runs from 1 to D sn (Schmidt number). An operational 
interpretation of eq. Hoi is as follows: measurement 011 E 
by the basis set {£, E } would result in S assuming state 
with probability A^. Now if <S, is the entropy of 
then, we can write directly from eq. [2J 

* — ^i— 1 

This is the desired expression for the LECM. 



This localization procedure may appear to face some 
problems when p s has degenerate eigenstates. Most of 
the time this difficulty can easily be resolved by the use 
of conserved quantities and the symmetries of the system 

We now discuss the important cases when the LECM 
can be shown to be optimal. We note that when mea- 
surement is performed by eigenstates of RDM (p E ), 
kij = kji = and the eq. [14] reduces to 



Tr Ay log 2 p^ (Q = Tr A* log 2 p^ (£.). (17) 

For spin- 1/2 systems, with both the parts Si and 
S2 taken to be single-sites, use of conserved quantity 
(here Z-componcnt of total spin) and parity or C2 sym- 
metry (for finite systems; for translationally invariant 
systems this will be automatically satisfied) leads the 
four eigenstates of p s into the following form: | tt)j 
^=(| t4) ± I 4t» and I ii), with t (i) being up (down) 
spin. Any two of the 2x2 matrices representing these 
states are seen to satisfy the condition given in eq. [T7] 
Since entropies of first and fourth states are zero, and 
second and third states (a singlet and a triplet) are one, 
the LECM in this case simply becomes, S = A s + A t , with 
A s (A t ) is the eigenvalue of p s corresponding to singlet 
(triplet). Note, after measurement, two sites will be in 
the statistical mixture of all the four eigenstates of p s 
and if X s = A t , then we will not be able to extract any 
useful entangled pair from the ensemble. 

Numerical Result.- We study the LECM (S) between 
two sites of a frustrated antiferromagnctic Hciscnberg 
chain (j'j — j 2 spin- 1/2 model [IH). Two sites are placed 
symmetrically as in fig. [Tp; this arrangement makes S to 
be an optimal. 

Its behavior against the distance between the sites (R) 
for different values of ji is shown in fig. [5] 

We see that, in the Neel phase (j'2 < 0.5), S falls with 
increasing R and reaches a constant value at large R. 
With increase in R, all the four eigenvalues of the p s be- 
come equal (hence A s = At), which results in the LECM 
assuming a value of 0.5 (which we call residual value 
or S r ). In fact, this particular value of the LECM is 
obtained if we take two sites (symmetrically) one each 
from two totally separate chains (unentangled) and per- 
form canonical measurement on remaining parts of the 
chains. Physically this implies that, in case of a sin- 
gle chain (where sites are connected), when R is large, 
two sites become unentangled, i.e., we can not localize 
'extractable' or useful entanglement between them by a 
canonical measurement. Keeping this in mind, we there- 
fore, can quantify actual extractable entanglement in our 
localization process as AS = S — S r . Any positive value 
of the quantity AS will give us the actual 'gain' in our 
localization process. In case of our single chain, the quan- 
tity AS falls with increasing i?, which can be understood 
qualitatively by the Valence Bond (VB) theory [Til [l5|. 
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FIG. 2. For the groundstate of a ji — j2 spin chain with 24 
sites, the LECM (5) against distance between two sites (R) 
is shown for different j2 values. The special case j2 — is just 
result for an antiferromagnetic spin- 1/2 Heisenberg chain. In 
the inset, result is shown for > 0.5 (the spiral phase). 

A groundstate can be expressed by linear combination 
of many VB basis states where basis states with nearest 
neighbor bonds (a bond represents an entangled pair) 
contribute more towards the groundstate compared to 
the ones with distant neighbor bonds. This says why 
for large R two sites become decoherent or unentanglcd. 
This fact naturally leads us to the notion of entanglement 
length (^e), which is the typical length scale upto which 
it is possible to localize useful or extractable entangle- 
ment between two sites. If fall in AS with increasing R 
is assumed to be exponential in nature, we can define £e 
as [H], [TU, — » — ln ^ s for large R. Since we have a 
finite system (24 sites), we can use two values of R and 
corresponding values of AS to estimate the value of £e 
in the following way: & = - lnA f jl^f' 52 ■ We take R = 
7 and 11 for this purpose (this particular values are cho- 
sen to avoid odd-even effect of a finite chain [HI), and 
calculated values of £e as a function of j'2 can be seen 
from fig. [3] 

The fall in the value of £g with increasing value of ji 
is not unexpected, as the contribution of VB basis states 
with long bonds decreases with the increasing value of 
j2- This is supported by the fact that at the Majumdar- 
Ghosh (MG) point (j'2 = 0.5) groundstate has only near- 
est neighbor bonding. 

The large oscillations in the value of S for j2 > 0.5 
(fig. [2]) can be understood by the fact that the phase of 
the system in this range is spiral in nature. The degree of 
entanglement between two sites depends on the relative 
phase factor between the sites. This is why for some 
distances the value of S is very low. 

Conclusion. -hi this letter, we have derived a simple 
but generic optimality condition, which would be help- 
ful in finding optimal values of the entanglement local- 
ized between two disjoint parts of a many-body system 
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FIG. 3. The entanglement length £e is shown here as a func- 
tion of J2 • 



by doing measurement on the remaining part of the sys- 
tem. Besides, we also discussed how it can be useful 
in studying mixed states of a general bipartite system. 
We further discussed a canonical way of localizing entan- 
glement which in some important and not-too-restricted 
cases (shown for spin- 1/2 systems) gives optimal value 
of the localized entanglement. The entanglement lo- 
calized by canonical measurement (or LECM) is oper- 
ationally meaningful and easy to calculate in practice. 
Unlike other measurement-based quantifications of en- 
tanglement, the LECM docs not require any demanding 
optimization procedure. Another important advantage of 
this LECM is that, since it does not depend on arbitrary 
choice of BSM, it provides a general framework for com- 
parative study of different types of quantum many-body 
systems. We studied a j\ — ]2 spin model to demon- 
strate the behavior of LECM. In this context, we also 
discussed extractable or useful part of LECM and de- 
fined an entanglement length scale upto which one can 
localize extractable entanglement between two sites. It 
may be stressed here that, all the concepts in this let- 
ter are quite general, virtually applicable to any kind of 
quantum many-body systems. 
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